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NUMERICAL INTEGRATION ORBITS AND BROUWER 


AND MODIFIED BROUWER ORBITS 


By 


Hans G. Hertz 


abstract 


An orbit constructed by means of Brouwer's theory (AJ 64.378.1959) was 
compared with a numerical integration orbit started with the same initial con- 
ditions and values for the coefficients of the earth's potential. 


Two sets of comparisons were made: one every 2 minutes for a 10-day arc, 
and one every 4 hours for a 60-day arc. 


The differences 'Brouwer- Integration 1 were appreciable but could be greatly 
reduced by modifying the computation of the short-period terms and by applying 
second order short-period terms to Brouwer’s expression for the semi-major 
axis as derived from Kozai's theory (AJ 67.446.1962). A further reduction was 
obtained by basing the mean motion of the mean anomaly on Kozai. 

The differences 'Modified Brouwer-Integration’ were represented as sums 
of polynomials of degree zero or one and Fourier series. These analyses re- 
sulted in further corrections to Brouwer. New comparisons after the applica- 
tion of these corrections which are empirical were made. 
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NUMERICAL INTEGRATION ORBITS AND BROUWER 
AND MODIFIED BROUWER ORBITS 


1. INTRODUCTION 

There are two fundamentally different methods of arriving at the quantitative 
description of the motion of an artificial satellite, i.e., at finding the position of 
the satellite at any given time. 

One approach is by using the method of numerical integration whereas an 
alternative is a representation of the motion by an analytical theory. The ad- 
vantages and disadvantages of these two alternatives are well known and need 
not be discussed here. But I feel there has been less discussion of the relative 
accuracy of these two methods. 

Also to the answer of this question two approaches may be considered. 

A set of observations could be represented by an orbit obtained by a numeri- 
cal integration with subsequent differential correction and, alternately, by an 
orbit derived from an analytical theory with an adjustment of the parameters of 
the theory to the observations included. The accuracy of the representation of 
the observations by the two orbits could then be compared. Because of errors 
of observations some theoretical shortcomings in either orbit may not always 
be easily recognizable. 

This criticism would not apply to the second approach which is the one 
selected here and which has, however, shortcomings of its own to be mentioned 
below. 

In view of the ’highest accuracy' of the numerical integration since it is 
subject to fewer approximations and inaccuracies than a general theory when 
carried out with a suitable interval a numerical integration orbit constructed 
with a certain set of initial conditions and a set of values for the coefficients of 
the earth’s potential will be taken as a ’standard orbit.’ An orbit will then be 
constructed with the help of a specific analytical theory, namely Brouwer's 
theory (Brouwer 1959) corresponding to the same initial conditions and the same 
values of the earth's potential as used in the numerical integration orbit. The 
deviation of the Brouwer orbit from the numerical integration orbit, the standard 
orbit, then gives a measure of the accuracy of the Brouwer orbit with the numer- 
ical integration as standard. The disadvantage of this method is that it is not 
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obvious whether any discrepancies which might show up are due to errors in the 
Brouwer theory or in the numerical integration or in both. Additional tests will 
be necessary to reduce or increase the probability for one or the other possibility. 


it will be seen that the results of the comparison may be varied by using 
modifications of the Brouwer theory. 


2. MATHEMATICAL MODELS 


As indicated in the introduction two mathematical models to represent satel- 
lite motion will be used in this investigation. 


The first model is a numerical integration. It is carried out by means of a 
computer program using the Cowell scheme. The nucleus of the program has 
been designed by Juergensmeyer (19621) but I have made many peripheral 
modifications and additions which, however, do not affect the mathematical 
model. 

The second model is based on the Brouwer theory (Brouwer 1959) which 
computes secular terms to the second order and periodic terms to the first 
order. A program to compute satellite positions according to Brouwer has been 
assembled by Rep ass and Chaplick (1964). I have again modified this program 
(Hertz 1969). The modified program allows the computation of satellite positions 
according to Brouwer as published by him with optional modifications and options 
to be included or excluded. Details as to the modifications and options have been 
described (Hertz 1969) but will be briefly stated here. 

The long-period terms can be omitted or included. The same applies to the 
short-period terms. If they are to be included they may be either computed from 
the secular portions e" ,1” as suggested by Brouwer or they may be computed from 
the long-period portions e' ,1' as would be more accurate on theoretical grounds. 
The lowest order terms in a neglected by Brouwer are the second-order short- 
period terms there being no long-period terms in a. Those short^period terms 
may affect the mean motion of the mean anomaly and thus the mean anomaly it- 
self. The program has an option to include them. They are based on expressions 
given by Kozai (1962) in his second order theory. If desired, the mean motion of 
the mean anomaly may be computed according to Brouwer or a reduction may be 
added which would bring it closer to the value based on the Kozai theory. Finally, 
specified Fourier series may be added to the osculating values of the Keplerian 
elements a, e, I,fi , <a, M at each data point before being used to compute a position- 
velocity vector for the data point. These Fourier series will be referred to as 
supplementary perturbations. 









Other options not being used need not be discussed here. The expression 
for the final value at a time t of an element a for use in computing the satel- 
lite’s position and velocity is of the form 

a = (T„ + 8 cr + 8 cr 

0 p s 

where cr Q is the secular portion, 8 p cr the periodic portion, and S 3 °- the value of 
the sum of the optional supplementary perturbations. In case of a ~ a, e, or I 
the secular portions are constants and in the case of cr - fi, co, or M they are 
linear functions of the time with constant coefficients. 


3. ORBITS CONSIDERED IN THIS INVESTIGATION 


For the purposes of this investigation we understand an orbit to be a set of 
orbital data generated in a computer run for a certain set of data points constructed 
from a specified set of coefficients of the earth’s potential, a specified mathe- 
matical model, and a specified position-velocity at a specified epoch. 

The orbital data considered here for each data point are 12 quantities, viz 
the rectangular equatorial coordinates and velocity components, and the osculating 
elements at the data point: 


a : 
e = 

I = 

n = 

Oil - 

M = 

Two orbits constructed with the same set of coefficients of the earth's 
potential, the same mathematical model, and the same position-velocity vector 

for the same epoch but with different sets of data points are considered to be 
different orbits but are, of course, expected to have the same values for the 
orbital data for any data point in common. 

For all orbits considered here we use with one exception one of two sets of 
data points: one consists of all 2 m points from the epoch to 10 d from the epoch, 
i.e. 7201 points; the other consists of all 4 h points, from the epoch to 60 d from 
the epoch, i.e. 361 points. In one case to be pointed out below the set of data 
points consists of all 40 m points from the epoch to 60 d from the epoch, i.e. 

2161 points. 


semi-major axis 

eccentricity 

inclination 

right ascension of the ascending node 
argument of perigee 
mean anomaly 


with respect to 
equator and equinox 
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Since we are interested in different mathematical models and two sets of 
data points we must consider several orbits between which comparisons are to 
be made. Altogether we consider eighteen orbits of which six are numerical 
integration orbits and twelve Brouwer or modified Brouwer orbits. They are 
described in Tables 4 and 5 and the text below. 

The same coefficients GM, J , J , J , J of the earth's potential (see Hertz 
1969, p. 1) 3 

U =-£ 


namely those given in Table 1 are used for all orbits except for orbits 5, 6, 17, 
and 18. However, also in these orbits the same value of GM as in the other 
fourteen orbits is used. The four orbits 5,6, 17, and 18 have been constructed 
for a special purpose discussed in section 8 which requires the use of different 
sets of values for J 3 and . 

All eighteen orbits use the same epoch and are constructed with the same 
position-velocity vector defined by the data in Table 2. Since in addition the 
same value of GM is used for all eighteen orbits the values oi the osculating 
elements a, e, I, 0, a, M at the epoch are the same for all the orbits. The 
values are given in Table 3. 


£ j„ (f 

n=l 


Jt =o, 


The fourteen orbits for which all coefficients in Table 1 are used consist 
of seven pairs, viz. orbits 1,2; 3,4; 7,8; 9,10; 11,12; 13,14; 15,16, Both orbits 
of each pair except in the case of orbits 1,2, and 3,4 have been constructed with 
the same mathematical model but one orbit uses the 2 m data points and runs for 
10 d whereas the other one uses the 4 h points and runs for 60 d . Also orbits 
1,2, and 3,4 have been constructed with the same mathematical model but the 
2 m and 40 m data points were used. Orbits 1 and 3 run for 10" and orbits 2 and 4 
for 60 d . 

Table 4 describes the numerical integration orbits. 

) 

Orbits 1 and 2 use the same mathematical model since they are both l m 
integrations. They will be used as the ’standard orbits’ for the 2 m and 4 h points 
and should, of course, be in complete agreement for any data points common to the 
2" and 4 h set. 

Orbits 3 and 4 being 0 n V5 integrations use the same mathematical model 
and will be used to test the reasonableness of the assumption to consider orbits 
1 and 2 as standard orbits. 
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Orbits 5 and 6 are used for a special investigation in section 8. 

From the coordinates and velocity components obtained in all numerical inte- 
gration orbits values for the osculating elements a, e, I,n , co , M were obtained. 

Table 5 describes the Brouwer and modified Brouwer orbits used. They 

. consist of the five pairs 7,8; 9,10; 11,12; 13,14; 15,16 to which the statements 

made above apply, and orbits 17 and 18 used in the discussion of section 8, 

. For orbits 7 and 8 the short-period terms have been computed with e n ,I” as 

suggested by Brouwer. No second-order short-period terms in a, no reduction 
to Kozai's mean motion, and no supplementary perturbations have been included. 
In view of what was said concerning this mode (Hertz 1969, p. 8) it would not 
be realistic to consider orbits with the short -period perturbations computed 
with e n ,I ,f , with second-order short-period terms in a or the reduction to 
Kozai T s mean motion included. 

In orbits 9-18 the short-period terms have been computed with e r ,I f which 
is the correct way on theoretical grounds (Hertz 1969, pp. 7, 8). 

In orbits 9-10 no second-order short period terms in a have been included, 
in orbits 11-18 they are present. The reduction to Kozai’s mean motion has been 
omitted in orbits 9-12 and 17-18 and included in orbits 13-16. The amount is 
+0?000 000 268 727 per hour for orbits 13 and 14 and +0?000 000 268 730 per 
hour for orbits 15 and 16. Supplementary perturbations have been included only 
in orbits 15,16 which will be discussed in section 7. 

Since the values of the osculating elements at the epoch are to be the values 
given in Table 3 it is these values which determine for each orbit and each a the 
value of cr Q at the epoch, i.e. the integration constant corresponding to cr, i.e. 
the mean value of cr , Its value is given by 

cr - cr - (S cr) -(Scr) . 

0 osc ' p ''epoch v s 'epoch 

where cr is the appropriate value from Table 3. and where (S^cr) and 

osc ' P 'epoch 

(S s cr ) e . poch are the values of S p o- and S s cr at the epoch. The value of cr Q there- 
fore changes with the mathematical model used. It therefore will have different 

• values according to whether the short-period terms are computed with e",I n or 
e*. I f or whether the 2 nd order short-period terms in a are included or not. The 
inclusion or non- inclusion of the reduction to the Kozai mean motion, on the 

* other hand, has no effect on cr Q . if the change from one mathematical model to 
another causes the change of do- in (S p ^) epoch + (S s °-) epoch then the value of cr 
will be changed by -dcr. 
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4. ORBIT COMPARISONS 


An orbit comparison is understood in this investigation to be the formation 
of differences of the form 


Aq = q 2 (t) - qj(t) 

where q x and q 2 are the values of the same orbital quantity q in orbits O i and 
0 2 respectively for the same time or data point t. 


These differences are formed for all data points common to both orbits. 
They are formed for the twelve quantities listed on p. 3 and are denoted by 


Ax, Ay, Az, Ax, Ay, Az, Aa, Ae, AI, AQ, Aa>, AM, 

and are said to have been formed in the sense 'Orbit 0 2 minus Orbit Oj , f They 
are functions of the time. 


These orbit comparisons will be made for various pairs among the eighteen 
orbits. The results of some of the comparisons will be used for further analyses 
to be discussed in section 7. 

For each orbit comparison made we will tabulate in Tables 6-9 twelve 
values one corresponding to each of the twelve quantities involved in the orbit 
comparison. If q a (t) and q 2 (t) are the values of q, where q is one of the 
twelve quantities, at time t in orbits and 0 2 respectively then we list the 
maximum value of | q 2 (t) - q x (t) | as t assumes all data points common to 
the two orbits. This value will be briefly called the maximum absolute value 
of the differences of q in the comparison of orbits O x and 0 2 and will also 
be denoted by Aq. 

Tables 6-9 contain the maximum absolute values of the differences of the 
twelve orbital quantities in all orbit comparisons. They will be discussed in 
sections 6 and 8 t 


5. THE NUMERICAL INTEGRATION ORBITS 

The numerical integration orbits 1 and 2 are the standard orbits for the 
2 m and 4 h data points. Since they have been constructed with the same mathe- 
matical model they should be in complete agreement for all data points in com- 
mon for the two orbits. Thus orbits 1 and 2 should agree for all 4 h points in the 


first 10 days. The comparison, actually made for every 40 m , showed complete 
agreement except in the inclination where the largest discrepancy is about 10 ~ 12 

degrees. 

No analytical determination of probable or maximum errors of the coordi- 
nates at the data points in the two standard orbits 1 and 2 has been made. In order 
to obtain some idea of the accuracy of the coordinates and the other nine orbital 
quantities in these orbits they have been compared with orbits 3 and 4 which are 
orbits with a 0T5 integration interval. The results are given in Table 6 and will 
be partially used in section 6. Graphs are presented in Figures la-f and 2a-f. 


6. COMPARISON OF THE BROUWER AND MODIFIED BROUWER 
ORBITS WITH THE NUMERICAL INTEGRATION ORBITS 


The results of the comparisons of the Brouwer and modified Brouwer orbits 
with the standard numerical integration orbits 1 and 2 are shown in Tables 7 and 
8 and Figures 3-12. The results in Table 7 and Figures 3-7 concern the 2 m data 
points while those in Table 8 and Figures 8-12 refer to the 4 h data points. 

The deviations of the unmodified Brouwer orbits 7 and 8 from the respective 
standard numerical integration orbits 1 and 2 are appreciable. They may be due 
to any one of the following four causes: 

1. an appreciable truncation error in the nu me rical integration program 
used in this investigation 

2. other types of errors in the numerical integration program 

3. an appreciable truncation error in the Brouwer theory 

4. errors in the compu. c program implementing the Brouwer theory and 
used in this investigation. 

Since the values given in Tables 7 and 8 for orbits 7 and 8 are very much 
larger than the maximum values of the differences between the 1“0 and 0T5 numer- 
ical integration orbits 1 and 3 and 2 and 4 given in Table 6 the first two of the 
possible causes for the discrepancies appear to be less likely than the other two 
although they are not completely ruled out. No effort will be presented in this 
paper to determine the actual cause of the discrepancies. However, in section 7 
we will consider the implications of the third cause. 

It is very puzzling that the maximum deviations from the standard orbits 1 
and 2 become larger as the computation of the short-period terms with e",!" is 
replaced by that with e', I f although the latter mode appears to be more prefer- 
able on theoretical grounds (Hertz, 1969, pp. 7, 8). No explanation can be offered 
here. 
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However, if the computation of the short-period terms with e’, I’ is adopted 
then the introduction of the second-order short-period terms in a greatly reduces 
the discrepancies, i.e. deviations from the standard orbits 1 and 2 in case of most 
of the twelve orbital variables. A further reduction is observed if the reduction 
to Kozai's value for the mean motion of the mean anomaly is included, at least if 60 
days are considered. Finally, an additional reduction except in the case of the ec- 
centricity for the 60 d run is achieved if the empirical orbits 15 and 16 to be dis- 
cussed in section 7 are used. The slight increase in case of the eccentricity does 
not seem to be significant. 

7. ANALYSIS OF THE DEVIATIONS OF ORBIT 14 FROM ORBIT 2 


Although we are not deciding which of the four causes stated on p. 7 are 
actually responsible for the deviations of orbit 8 from orbit 2 we shall investi- 
gate what the implication would be if the actual cause is No. 3 listed. 

The differences between these two orbits and also those between orbits 14 
and 2 would then be expected to be Fourier series in case of a, e, I and the sum 
of a linear function of the time and a Fourier series in case of Q , or, and M. The 
arguments in these Fourier series would be linear combinations with integral 
coefficients of M and a> where M and co are the secular portions of M and co as 
obtained in the modified Brouwer orbit 14. If such Fourier series and sums of 
linear functions of the time and Fourier £ eries would be applied to orbit 14 with 
the proper signs we would expect to obtain an orbit which would be much closer to 
orbit 2 than orbit 14 and orbit 8 are. 


An attempt was made to find the stated representations for the differences 
between orbits 2 and 14. For this purpose a computer program which cannot be 
described in detail here has been used. This program has the following capabilities: 

The values of a function of a single variable t given for a set of equidistant 
values of t are represented by the sum of a polynomial in t of specified degree 
and a Fourier series. The arguments are linear functions of t. The program 
determines the number of arguments to be used, the coefficient V of t for each 
argument and the coefficients of the terms cos \ t and sin X.t the presence of 
which are implied by the selection of the value k by the program. The k's 
will be called the free frequencies in the analysis. If desired the function can 
be assumed also to have Fourier terms with specified or prescribed frequencies. 
The program will determine the coefficients of the associated cosine and sine 
terms. 


For the investigation on hand there would be required six analyses: one 
each for the differences of a, e, I, ft , co, M in orbits 2 and 14 for every 4 h from 
the epoch to 60 d after the epoch. Ideally, the frequencies of all arguments of the form 
iM + jcd where i, j assume all combinations of integers which could reasonably be 
expected, would be taken as prescribed frequencies. However, the program allows 
the use of only twelve prescribed frequencies, which would not always provide 
enough frequencies since it was not a priori known which frequencies would actually 
occur. The difficulty was circumvented in the following way: in a first iteration 
only the frequencies co , 2 co, . . . , 6di were prescribed since they were expected to 
belong to important long-period arguments on theoretical grounds . If k is any free 
frequency found then a term with cos k t or sinK-t has at the points t = in At* the 
same values as the term cos(3t or sin/3t with the same coefficients and with 
,6 = k + 2-7rk/At, k being an integer. In other words, the actual frequency of the 
term could be/3 = k + 27rk/A t instead of k with some integral value of k which 
cannot be determined from the data. But in all eases a reasonably small integer 
k could be found such that/3 is close to an integral multiple of M. Therefore, since 
co is small compared with M two integers i and j can be found such that 

j/3 - (i M + j co) | < i-jwf 

The analysis now was repeated with introducing prescribed frequencies iM + j co 
corresponding to each free frequency. Because of differences up to (1/2) | ^between 
the free and corresponding prescribed frequencies slightly different coefficients 
might result in the new analysis. The new analysis might generate new free fre- 
quencies. The procedure was repeated until no more free frequencies arose. 

Seldom more than twelve prescribed frequencies were required. If it happened 
it was possible to drop one or two frequencies jco with small coefficients to 
accommodate the extra frequencies. In view of the higher level iteration proce- 
dure described below this is not expected to introduce any serious error. 

The results of such analyses led to Fourier series or sums of linear functions 
of the time and Fourier series which, when applied with the proper signs to a, e, I, 
n , co , M of orbit 14 would give an orbit which might be expected to be closer to 
orbit 2 than orbit 14 is. A slight modification must be mentioned. The constant 
terms obtained in the analyses were not used. The equation on p. 5 shows that 
changes within certain bounds in the constant terms of the supplementary per- 
turbations have no effect on the values of the supplementary perturbations for any 
t since such changes are accompanied by opposite changes in the values of the 
corresponding mean elements. Therefore the constant terms in the supplementary 
perturbations for ft, co , M have been put equal to 0. In case of a, e, I they may 
uniquely be determined in such a way that the changes in the motions of ft , co, M 
as required by the analyses may be realized. 

* m = 0, 1, ..... 
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With the new orbit thus obtained we repeat the cycle. Altogether four such 
cycles were passed through. The final orbit, orbit 16, was considered to be as 
close as possible to orbit 2. For reason of shortage of time no further cycles 
were carried out but it is believed that orbits resulting from additional cycles 
would change the deviations from orbit 2 only insignificantly. Orbit 15 is then 
the orbit computed with the same position-velocity vector and the same supple- 
mentary perturbations for the 2 m data points. It was already stated on page 8 
that the deviations of orbits 15 and 16 from orbits 1 and 2 were smaller than 
those for orbits 13 and 14, and, of course, smaller than for orbits 7 and 8. 

Table 10 shows the supplementary perturbations which must be applied to 
both orbits 13 and 14 in order to obtain orbits 15 and 16. To obtain orbit 16 
from orbit 8 we have to replace the computation of the short-period terms with 
e",I" by that with e’,I\ add the second-order short-period terms in a, the re- 
duction to Kozai's mean motion, and the supplementary perturbations of Table 10. 
The total of these reductions then would represent the truncation error of the 
Brouwer theory if it were actually the cause of the deviations of the unmodified 
Brouwer orbit 8 from the numerical integration orbit 2. 

The hypothesis could be tested if a more accurate theory such as Kozai's 
theory {Kozai 1962) would be compared with a numerical integration orbit as 
has been done here. 


We shall conclude this section with two remarks. 

The supplementary perturbations in Table 10 for a contain terms with 
arguments jco, i.e. long-period terms. Such terms are not present in Brouwer's 
and Kozai's theories. It has not been determined whether there could be higher 
order terms of this type. Possibly the analyses should have been repeated with- 
out using co, 2 co, . . . 6 co as prescribed frequencies. 

Orbit 16 is closer to the integration than orbit 8. Orbit 14 is also closer 
but not as close as orbit 16. But orbit 14 can be constructed by 'analytical' 
means whereas orbit 16 is an empirical orbit and moreover requires the 
knowledge of the numerical integration orbit. The advantages of analytical 
orbits would be defeated if one were to try to improve on unmodified Brouwer 
orbits by applying all corrections which in our case have led to orbit 16. If we 
want to adhere strictly to analytical processes we can short of introducing a 
more accurate theory improve the unmodified Brouwer theory by: 

1. Replace the computation of the short-period terms with e",I" by the 
computation with e' 

2. add the second-order short-period terms in "a" 

3. add the reduction to Kozai's mean motion of the mean anomaly. 
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It is hoped to present a more through interpretation of the orbit comparisons 
in a future report. 


8. COMPARISON OF ORBITS WITH SOME ZERO HARMONICS 


The numerical integration orbits 5 and 6 are identical with numerical inte- 
gration orbit 2 except that in orbit 5 J 4 and in orbit 6 J 3 and J 4 have been put 
equal to zero. The modified Brouwer orbits 17 and 18 are related similarly to 
modified Brouwer orbit 12. The results of the comparison of orbits 17 and 18 
with orbits 5 and 6 are shown in Table 9. These values should be compared with 
the results of the comparison of orbits 12 and 2 shown in Table 8. The elimination 
of non-zero values of J 4 and J 3 brought some reductions in the differences 
' Br ouwer-Integr ation' but also some increases. It is thus hard to say whether 
the discrepancies for orbit 12 are due to J 2 or the higher harmonics. 
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Table 1 

Coefficients of the Earth’s Potential and Equatorial 
Radius of the Earth Used in Orbits 1-18* 

GM 3.98603200 * 10* km* seconds- 2 

J 2 +1.0823 • 10 ' 3 

J 3 -2.3 • 10" 6 

J 4 -1.8 *10~ 6 

Jg 0. 0 

R 6378.165 km 


*in orbits 5, 6, 17, 18 J 4 = 0 

in orbits 6, 18 J 3 = 0 


Table 2 

Epoch and Position-Velocity Vectors 
Used for all Orbits 

Epoch 1966 November 29 0 h 00 m 00 s ET 

Position-velocity vector with respect 
to the equator system 

+4 529 467.000000 "1 

+5 724 243.000000 > meters 

-3 803 180.000000 J 


-18 438 869.000000') 

+14 616 528.000000 > meters per hour 
- 8 988 297.000000 J 
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Table 3 

Osculating Elements at Epoch 
for all Orbits 


Epoch 1966 November 29 0 h 00 m 00 s ET 

a 8 321 016.341040682 meters 
e 0.1659571587412389 

I 32f8537259317278 

n 17 7°8 6 18 8 9 142 1345 

co 142?5592123598749 

M 76^7687107985003 


In the actual runs there may be values different from those given here by 
a few units of the last place given. 


Table 4 

Numerical Integration Orbits 


Epoch and position-velocity vector at epoch from Table 2 


Ser. No. 

Coefficients of the 
Earth’s Potential 

Integration 

Interval 

Total 

Length 

Ephemeris 

Interval 

1 

Table 1 

1*0 

10 d 

2 m 

2 

Table 1 

1.0 

60 

40 

3 

Table 1 

0.5 

10 

2 

4 

Table 1 

0.5 

60 

40 

5 

GM, J 2 , J 3 from Table 1 
J 4 = J 5 = 0 

1.0 

60 

240 

6 

GM, J from Table 1 

T T T n 

1.0 

60 

240 


3 4 5 



Table 5 

Brouwer and Modified Brouwer Orbits 
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Table 6 

Maximum Absolute Values of Differences in Comparison of the 
numerical integration orbits (orbits 3 and 4) with the 1T0 
numerical integration orbits (orbits 1 and 2) 


Orbit 3 - Orbit 1 Orbit 4 - Orbit 2 


Ax ^ 0.09 2.88 

Ay > in meters 0.09 2.49 

Az J 0.06 1.45 
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Ay W in meters per hour 0.3 9.4 

Az J 0.2 6.0 
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A Cl) 


>• 
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2 m points 4Q m points 
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Table 7 

Maximum Absolute Values of Differences in Comparison 
of Brouwer and Modified Brouwer Orbits with l m Numeri- 
cal Integration Orbit using 2 m Points 
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Table 9 

Maximum Absolute Values of Differences in Comparison 
of Modified Brouwer Orbits and lUO Numerical Integra- 
tion Orbits with Some Zero Harmonics Using 4 h points. 

Orbit m - Orbit n 

Orbits m are modified Brouwer orbits with the short period terms com- 
puted with e\I t and second order short-period terms in a included. No 
reduction to Kozai's mean motion and no supplementary perturbations 
were included. 


J 2 is taken from Table 1 
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Ax •'j 
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A z J 

I 

150.8 

109.3 

Ax 
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A y 

► in meters per hour 
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391 

Az J 
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370 

Aa 

in meters 

0.3 

0.2 

Ae 

in units of 10 “ 7 

25 

16 

AI 1 


10 

4 

AQ 1 

> in units of 0?00001 

63 

46 

tso) 

139 

128 


Table 10 

Supplementary Perturbations to be Applied to Orbits 13 and 14 
Arguments are 0 at JD 2439488.5 
Motions of arguments are of the form iM + j <i>, where M and & 
are secular motions of M andeo in orbit 13 or 14 respectively 
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Figure Ic-Plot of Data from Comparison of Orbit 3 Minus Orbit 1 Inclination. 
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Figure 3a-Plot of Data from Comparison of Orbit 7 Minus Orbit 1 Semi Major Axis. 
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Figure 3c-Pfot of Data from Comparison of Orbit 7 Minus Orbit 1 Inclination 
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Figure 3e-Plot of Data from Comparison of Orbit 7 Minus Orbit 1 Argument of Perigee. 
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THE ABSCISSAE ARE MINUTES FROM THE EPOCH 

Figure 4c^Plot of Data from Comparison of Orbit 9 Minus Orbit 1 Inclination. 






Figure 4e-Plot of Data from Comparison! of Orbit 9 Minus Orbit 1 Argument of Perig 
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Figure 5c-Plot of Data from Comparison of Orbit 11 Minus Orbit 1 Inclination 





Figure 5e-Plot of Data from Comparison of Orbit 11 Minus Orbit 1 Argument of Perigee. 



Figure 5f-Plot of Data from Comparison of Orbit 11 Minus Orbit 1 Mean Anomaly. 
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Figure 6a-Plbt of Data from Comparison of Orbit 13 Minus Orbit 1 Semi Major Axis. 




Figum 6c-Plot of Doto from Comparison of Orbit 13 Minus Orbit 1 inclination. 




Figure 6e-Plot of Data from Comparison of Orbit 13 Minus Orbit 1 Argument of Perigee. 
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Figure 7a -Plot of Data from Comparison of Orbit 15 Minus Orbit 1 Semi Major Axis. 
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Figure 7d-Plot of Data from Comparison of Orbit 15 Minus Orbit 1 Right Ascension of Ascending Node. 
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Figure 7e-Plot of Data from Comparison of Orbit 15 Minus Orbit 1 Argument of Perigee- 








Figure 8b-Plot of Data from Comparison of Orbit 8 Minus Orbit 2 Eccentricity. 
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Figure 9a-Plot of Data from Comparison of Orbit 10 Minus Orbit 2 Semi Major Axis. 
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Figure ?c -Plot of Data from Comparison of Orbit 10 Minus Orbit 2 Inclination. 




Figure 9e-Plot ov Data from Comparison of Orbit 10 Minus Orbit 2 Argument of Perigee. 
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Figure 10o- Plot of Data from Comparison of Orbit 12 Minus Orbit 2 Semi Major Axis. 
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Figure ITa-Plot of Data from Comparison of Orbit 14 Minus Orbit 2 Semi Major Axis. 


REPRODUCIBILITY OF THE ORIGINAL PAGE IS POOR. 
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Figure 12c-Plot of Data from Comparison of Orbit 16 Minus Orbit 2 Inclination 





